Squeezed noise due to two-level system defects in superconducting resonator circuits 
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Motivated by recent surprising experimental results for the noise output of superconducting microfabricated 
resonators used in quantum computing applications and astronomy, we develop a fully quantum theoretical 
model to describe quantum dynamics of these circuits. Building on theoretical techniques from quantum optics, 
we calculate the noise in the output voltage due to two-level system (TLS) defects. The theory predicts squeezing 
for the noise in the amplitude quadrature with respect to the input noise, which qualitatively reproduces the noise 
ellipse observed in experiment. We show that noise enhancement along the phase direction persists for pump 
frequencies away from resonance. Our results also suggest that intrinsic TLS fluctuations must be incorporated 
in the model in order to describe the experimentally observed dependence of the phase noise on input power. 

PACS numbers: 03.67.-a, 03.65.Yz, 11 22. Gm, 85.25.-j 



Amorphous dielectrics contain weakly-coupled two-level 
system (TLS) defects QIS that are known to modify the com- 
plex permittivity at low temperatures, including the imaginary 
part which is responsible for loss |3 1. In Josephson qubits the 
purposeful amorphous dielectrics were found to contain TLSs 
which leads to decoherence |4 |. Similarly, the loss in super- 
conducting resonators is generally attributed to TLS defects 
which can arise from the native oxides on the superconduct- 
ing leads [51 or deposited amorphous dielectrics |E1[3. Su- 
perconducting resonators are used in applications of quantum 
computing [8 1 and as single photon detectors for astronomy, 
where noise can limit the detection performance 1 9 1 and may 
be attributable to TLSs |10l. Surprisingly, recent studies have 
revealed that the amplitude noise quadrature of resonators is 
limited by the experimental sensitivity near the vacuum noise 
limit ifTTIl . In contrast, noise in the phase quadrature is rela- 
tively large and decreases as the square root of the measure- 
ment power |[T2ll . A semi-empirical model was used to de- 
scribe this phase noise saturation, but a quantitative model of 
phase-amplitude noise asymmetry is still lacking 1 13 |. 

Previous related noise studies have considered the effects of 
TLS noise on superconducting qubits |T4l, and noise squeez- 
ing in various nonlinear systems including Josephson junc- 
tion parametric amplifiers (151 [El. micro-cantilevers ifTTl , 
and nanomechanical resonators coupled to a Cooper-pair box 
ifTSl . Atoms have also been used for four- wave mixing exper- 
iments to squeeze optical cavity states |[T9l . 

In this letter, we develop a theory for noise due to TLS 
defects in a superconducting resonator circuit. We compute 
the noise in the transmitted voltage for the system shown in 
Fig. [TJa) assuming that the lone source of stochastic fluctua- 
tions are the quadrature-independent white noise in the inci- 
dent voltages on the transmission lines. The theory predicts 
that squeezing in the transmitted noise generally occurs even 
if the dynamics of the internal resonator, including the TLSs, 
is fully deterministic. The results imply that intrinsic stochas- 
tic dynamics of TLSs is not necessary in explaining the ob- 
served squeezing in the transmitted voltage noise. Our work 
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FIG. 1: (a) A circuit diagram of the considered setup, (b) Actual 
image (taken from Ref. 6 ) of a notch-type aluminum LC resonator 
which can be modelled using (a). The components are labelled to 
clarify its correspondence to circuit in (a). 



shows a similarity between noise squeezing from TLS defects 
in these devices Il2l and squeezing of coherent light due 
to two-level atoms in quantum optics |20|. We also compute 
noise power for ofl'-resonant pump frequencies and confirm 
that enhancement in the phase noise and squeezing in the am- 
plitude noise persists away from resonance. 

We have two identical semi-infinite ID transmission lines, 
labeled - and +, which are attached to each other end-to-end 
at X = (see Fig. [TJa)). They are coupled to the internal 
superconducting LC resonator at x = through the coupling 
capacitor Q. The capacitor and inductor in the LC circuit 
have capacitance C and inductance L. Here, we take Q <^ C, 
and consider the case where TLS defects reside only inside the 
dielectric of capacitor C. A microwave source (not shown) 
sends in a pump signal from x = -oo, and we focus on the 
transmitted voltage noise on transmission line +. 

The Hamiltonian for the full system has three terms, H = 
Ho + ^TLS + Hji^s-R- Ho = ,-} + Hr modcls the two 



2 



transmission lines and the internal resonator, 



dx 0(mx) 



21 



2c 



^ 2L 2Cr 2C 



(1) 
(2) 



where m = {+, -} labels the left (-) and right (+) transmis- 
sion lines, and 0(x) is the unit step function. Qm(x) is the 
operator for total charge residing to the right of x on transmis- 
sion line m, Ql is the operator for total charge that has flowed 
through the inductor, and Qo = (Q- - Q+)\x=o denotes the 
total charge operator on capacitor Q. Pm and Pi are the con- 
jugate momenta for and Qi, respectively, and they satisfy 
[Qm(x),Pm'(x')] = ifidjnm'^ix - x') and {Ql.Pl\ = iK where 
m, = {+,-}. The transmission lines are modeled as conduc- 
tors with inductance per unit length capacitance to ground 
per unit length c, and characteristic impedance Z = 
1211 [22l . We assume we have N identical and independent 
TLSs, all with asymmetry energy and tunneling energy Aq. 
The TLS Hamiltonian is then given by T^tls = z + ^o^x. 
where Si = Tja=i ^ia is the collective spin operator which 
represents the N TLSs, and Sia = cr//2 with the usual Pauli 
matrices ct/. The components of the collective spin operator 
obey = ieijkSk- After diagonalization T^tls = ES^ 

where E = [A^ -h A^]^/^. Each TLS interacts with the uniform 
electric field inside the capacitor through its electric dipole 
moment p. We assume the dipoles fluctuate by making 180° 
flips between parallel and anti-parallel orientations with re- 
spect to the field. The field aff'ects the asymmetry energy 
and here we ignore the relatively small changes to the tun- 
nel barrier, similar to other treatments of TLSs derived from 
the tunneling model (IJ. In the diagonalized basis the inter- 
action between the TLSs and the field can then be written 
as Hji^s-R = gn[S, cot2^ + (S^ + S-)/2](Qo - QlX where 
the coupling constant g = 2|p| sin 2^/fiCd, d is the sepa- 
ration between the plates of the capacitor, ^ is defined via 
tan2^ = Ao/Aa, and ^- = ± iSy. We use |p| = 1 Debye, 
which is comparable to the dipole moment sizes observed for 
TLSs in amorphous Si02 at microwave frequencies and OH 
rotors in AlO;, ll23ll . 

For X 9^ 0, Qm(x, t) obey the massless scalar Klein-Gordon 
equation, and the solution can be written as a sum of right- 
{R) and left-propagating (L) components, i.e. Qm^xj) = 
QmR(x, t) + Qmdx, 0, where 
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V = {icY^I^ is the velocity of wave propagation and is the 
resonance frequency of the loaded LC circuit. The charge op- 
erators obey [^mrM^C'T'^^')] = i7lfll0.rZ)dnim'^TT'^{CD-(J^'). 

where r, = {R, L}. We also introduce the ladder operator for 
the inductor charge via Ql = qL+q\, where [qi, q^j} = fi/2ojoL 
and ojo = {LCT^'^. 



At X = 0, the quantum Max well-B loch equations for our 
problem read 

^AQ^R + Q^Dl^o^ (4) 

-(1-5.) 

LC^l-{Q^-Ql) = Cgn(S,coi2^ + (S^ +S-)/2)(6) 
[i{E/n) + T2]S- = ig{S,-coi2^S-)(Qo-QL) (7) 

(8) 
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Here, we have introduced phenomenological longitudinal 
and transverse TLS relaxation rates, Fi = and F2 = 
and Go = [Q-R + Q-L-Q^R-Q^ii ^? - 
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-(N/2) tSinh(E/2kBT) is the equilibrium expectation value for 
Sz when the TLSs are decoupled from the fields. Assuming 
we are in the low- temperature regime, i.e. «: E, we 

take ^ -N/2. The first three equations above are quantum 
generalizations of the Kirchhoff' voltage law, and the last two 
equations are quantum Bloch equations describing the cou- 
pling of the TLSs to the photon fields. 

We now perform a series of standard approximations which 
maps the above equations to a form often seen in quantum 
optics literature describing an ensemble of two-level atoms in 
an optical cavity 1201 . Since we have a high- 2 resonator we 
perform the Markov approximation, where we replace ^fco by 
^/D^ in the mode expansion Eq. ^ and extend the lower limit 
of the integral to negative infinity. We now go to a frame ro- 
tating at the pump frequency Q, and use Greek letters to de- 
fine operators in the rotating frame, i.e. 6^x0) = qmriOe^^^ 
e^{t) = qL{t)e'^\ t^it) = S^(t)e^'^\ and \{t) = S,(t). Per- 
forming the slowly-varying envelope approximation on the 
QlO) term and with the usual rotating- wave approximation, 
Eqs. Q-ljS]) can then be rewritten as 



zcM-mo^R-^L) 
-Lc(p} + 2i^id^eL 

{dt + /A + F2)£- 
25,£, + 2Fi(£,-sO) 



(9) 
§0 - LCc (p^ + 2i^dt) §l(10) 

%-eL + {Cgn/2)t- (11) 

igt,{eo-h) (12) 
ig\t-{el-el)-h.c.^, (13) 



where A = is the TLS detuning, and 6 [0,1] = 6-R+6-1- 

0+R + 0+L. We note that the voltage fields propagating along 
the transmission lines relate to the charge fields Omr through 

We solve Eqs. ([9])-(p3]) for A/^ » 1 where the dynam- 
ics of the quantum system can be described semi-classically 
|[24l . We treat each field (p(t) as a c-number and assume it 
has a steady-state mean part and a fluctuating part S(p(t) 
which models the noise arising from the incident voltage on 
the transmission lines. Assuming small fluctuations, we lin- 
earize Eqs. ([9|)-(p3]) with respect to these fluctuations about 
the steady-state value. 

The steady-state solution to the transmission amplitude 
v+r/v-r is plotted in Fig. [2ja) for real V-r (input voltage am- 
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FIG. 2: (a) Steady-state solution for the transmission with N = 100, 
and (i) = lOyuV; (ii) = lyuV; and (iii) v.r = 0.05//V. Here, 
AQ ^ 0.0003a;o. (b) Noise ellipses for v^r = 2yL/V, Q = Q„ and 
= 100 (red), 300 (green), 600 (blue), and 1000 (magenta). The 
dashed line is quadrature-independent noise with no TLSs. (c) Plot of 
population imbalance, 5 2, as a function of input power for N = 100. 



plitude) near the resonance frequency Q^- Blue points corre- 
spond to pump frequencies below and red points to those 
above Q^. The solutions are plotted for three different values 
of v-R, and the parameters used are Q = 0.0 IpF, C = 0.3pF, 
Z = 50Q, /o = oJo/27: = 6 GHz, Ti = 300ns, T2 = 30ns, 
d ^ 70nm and N = 100. We also assume the TLSs to be 
in resonance with the LC resonator and fix the TLS energy 
to E = fir throughout and take = Aq. We then find 
~ 0.98376i;o. In Fig. ^c), (TLS population imbal- 
ance) is plotted as a function of the input power, As P/„ 
increases approaches zero signifying TLS saturation. We 
see that saturation occurs for Pin > 10"^^W. Fig. ^sl) shows 
that the transmission deviates very little from the defect-free 
(no TLSs) limit for saturated TLSs. This is expected since 
only a small fraction of the total energy is stored in the TLSs 
when they are saturated. 

We now move on to obtaining the noise in the transmit- 
ted voltage. Solving the linearized equations the transmitted 
charge fluctuations can be written in terms of the two input 
fluctuations as 

A^K(aj)S©^R(aj) = A.K(aj)S©.R(aj) + A^^(a;)^0+L(^^), (14) 

where oj is the frequency deviation measured from the 
pump frequency Q, S&mA^) = i^OmT(ci>),SOl^^(-ctj)Y , and 
the elements of the coefficient matrices AmAoj) are given 
in the Supplementary Material. To study the quadrature- 
dependence of the noise in the transmitted voltage we intro- 
duce generalized voltage fluctuation variables Svmri^,^) = 
[e~^'^SvmT(oj)-\-e^'^Svl^^(-oj)]/2, where (f is the quadrature angle 
measured with respect to the real axis. We focus on the sym- 
metrized correlator of these fluctuations, i.e. S +r(oj, oj\(fi) = 
{^{Sv+r((jj,(p),SvIj^((jj\(P)}), where {A,B} = AB -\- BA. As- 
suming that the incoming fluctuations, Sv-R(t) and Sv+i(t), 
are characterized by a quadrature-independent white noise 
spectrum ^0 (for thermalized transmission lines it is s^^ = 




FIG. 3: Colour intensity plot of the normalized noise power 
s+R(aj,(p)/ So on the (oj, ^)-plane for pump frequencies on and away 
from resonance, V-r = 2fiV and N = 1000. Plots correspond to noise 
centred around (a) west; (b) northwest; and (c) north points, where 
north is defined as the top point on the resonance circle. 



27iflZflr coth(fiflr/2kBT)), we have 
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Using Eqs. (14) and (15) together with Sv+Rico^cp) ^ 
-2iZQr^6+R((jj, (f), S +r((jj, (jj\(p) = s+r((jj, (p)S((jj - oj^) can be 
straightforwardly obtained, (see Supplementary Material for 
details). In the following, we focus on the resulting spectral 
density of the noise, s+r((jij, (p). 

We first focus on the on-resonance case (where Im{v+/?} = 
0). Fig. [2jb) plots noise ellipses centred at this resonance 
point for v.r = 2/jY and various A^. Here, we are plotting the 
61; = component of the noise spectral density normalized by 
^0- The noise ellipse is defined such that every vector from the 
origin to a point on the ellipse makes an angle (f (quadrature 
angle) with respect to the positive real axis and has length 
s+r(<jj = 0,(fi)/so. In the defect-free limit, the noise ellipse 
is circular with radius 1 (the dashed line), which shows that 
in the absence of non-linearity the transmitted noise remains 
quadrature-independent. In the presence of TLSs, we obtain 
squeezing where the noise along amplitude (phase) quadrature 
is reduced below (enhanced above) the noise of the incoming 
fluctuations. Eccentricity of the noise ellipse increases as the 
number of TLSs is increased. 

Fig. |3] is a colour intensity plot of the normalized noise 
power s+R(oj,(f)/ So on the (cl), ^) -plane for on- and off- 
resonance pump frequencies. Here, V-r = 2/jy and = 1000 
are both fixed. We now use the convention where north cor- 
responds to the point at the top of the resonance circle in Fig. 
|2ja), and west to the left-most point on the circle and so on. 
Then, the three plots correspond to the noise centred around 
(a) west; (b) northwest; and (c) north. The results show that 
the major axis of the noise ellipse is always along the phase 
direction. Our results are consistent with experiments, where 
ffuctuations are also primarily observed in the direction tan- 
gent to the resonance circle lfT2ll . 

In Fig. |4ja) we plot the normalized excess phase noise 



power, ^exc(^) 



[s+R(ci),(p = 71 12) - so\Isq and nega- 



tive of the normalized excess amplitude noise, -s^^^icS) := 
-[s+R(aj,(f = 0) - so]/so, as a function of frequency devia- 
tion away from resonance, oj. Here we use v-r = 2/dV and 
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FIG. 4: (a) Frequency-dependence of the excess phase noise power 
and negative of the excess ampHtude noise power for = IjiN and 
= 1000. Roll-off at ^ cjjIcjJq - can be seen, (b) Normahzed 
excess phase noise at a; = and inverse of the internal quality factor 
as a function of input power (Q = Q^). 



can lead to observable noise fTSl . We reiterate, however, that 
squeezing phenomenon itself is present regardless of deter- 
ministic or stochastic nature of TLS dynamics. 

In conclusion, we have developed a theoretical framework 
applicable to circuits containing transmission lines, lumped 
circuit elements, and TLS defects, building upon previous 
work in quantum optics. It serves as a first step toward a 
quantitative theory for noise in these circuits. We found that 
quadrature-independent incident noise is generally squeezed 
once transmitted through a lumped circuit containing TLSs 
even when the latter evolve deterministically. Extensions of 
the model could allow for further quantitative results on noise 
due to TLSs, including the treatment of the standard tunneling 
model distribution of TLSs and incorporating intrinsic TLS 
fluctuations. 
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= 1000. Both noises roll off at the resonator bandwidth 
^roii/^o ~ 10""^- From the steady-state solution the quality 
factor of the resonator is estimated to be Qr ~ 9800, which is 
consistent with the plot. We see that the phase noise stays 
above the incoming noise while the amplitude noise re- 
mains below this value signifying squeezing. They both ap- 
proach ^0 as the frequency deviates sufficiently beyond the 
resonator bandwidth (see also Fig. |3]). For oj » 60)^11 we find 
^^xcM ~ and ~ {colco^)-'^-^. 

In Fig. we plot the loss tangent, QrJ for the pump 
frequency on resonance and = 100. This was found nu- 
merically by fitting the transmission amplitude to a circuit 
model as a function of [¥(1 + 2ext/Gint)]VZ ~ V^Z. Here, 
V is the voltage across capacitor C (as shown in Fig. [T]) 
and 2ext = 2(C + Cc)l{(ji>oZCl) is the external quality fac- 
tor. The above approximation holds because our resonator is 
over-coupled. At low amplitude the loss is constant, and at 
high amplitude the slope is ^ -1 as the spins are saturated. 
This is in agreement with the semi-classical theory for a sin- 
gle TLS type. For the standard tunneling model distribution 
the density of defects follows P(Aa, Aq) = Pq/Aq and the su- 
perposition of the loss from diff'erent TLSs give a slope of 
-1/2 in the high power regime 1 3 1, as observed in amorphous 
films. 

In the same figure the excess phase noise s^^ci^ = 0) is 
plotted. In the same high power regime as the loss tangent, 
the partially saturated TLSs exhibit phase noise with a slope 
of -1, which disagrees with the power dependence observed 
in experiment |9|. This suggests that a correct description for 
the phase noise power dependence necessitates the inclusion 
of intrinsic stochastic ffuctuations of the TLSs. Indeed, in- 
trinsic TLS noise causes dielectric constant fluctuations and 
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